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We argue that the cohinear factorization of the fragmentation functions in high 
energy hadron and nuclei cohisions breaks down at transverse momenta fcj- < Qs/ 9 
(yQ \ due to high parton densities in the colliding hadrons and/or nuclei. We calculate, at 

o 

CN) ' next-to-leading order in projectile parton density and to all orders in target parton 

(N : 

■ density, the double-inclusive cross section for production of a pair of gluons in the 
!>■ ■ 

\ scalar J^^ = 0^^ channel. Using the low energy theorems of QCD we find the 

Q^l inclusive cross section for vr-meson production. 
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;h ' Strong interactions at high energy are crucial 



y influenced by high density of gluons in 



113,3. 



A systematic description of parti- 



wave functions of colliding hadrons and/or nuclei 
cle production in high energy hadron reactions can be carried out in terms of quasi-classical 
solutions to the Yang-Mills equations on the light-cone Q] . The quasi-classical approxima- 
tion is valid for field modes which occupation number reaches the saturation limit oi ^ 1/ g. 
These modes correspond to gluons with transverse momentum less than the saturation mo- 
mentum Qg. The square of the saturation momentum is a measure of the two-dimensional 
color charge density in the hadron/nucleus wave function. As the collision energy increases, 
the quantum fluctuations increase the occupation number of modes with larger transverse 
momentum which manifests in increase of Qs- The corresponding high energy evolution of 
scattering amplitude is governed by the nonlinear Balitsky-Kovchegov evolution equation 



B.a. 



A comprehensive reviews of the physics of gluon saturation (Color Glass Condensate) 



can be found in Refs. 



71, |8j. A quasi-classical approach to inclusive processes of high energy 



QCD has driven a lot of attention due to its remarkable phenomenological success in small- a: 



DIS, p{d)A and AA collisions 



ifli. 



One of the challenges of QCD is to understand how the color degrees of freedom meta- 
morphose into hadrons. In the traditional perturbative QCD this problem is solved with 
the help of the coUinear factorization theorems which allow to separate the universal non- 
perturbative parton distribution functions, fragmentation functions, and the hard partonic 
sub-processes. The coUinear factorization theorems hold only if particle production is char- 
acterized by a momentum scale which is much larger than the typical momentum scale in 
hadron wave functions. At high energies, existence of a semi-hard scale Qs - which is an 
increasing function of energy and atomic mass A - does not allow application of factorization 
theorems at transverse momenta of the order of Qg or smaller ^. This is the kinematical 
region in which bulk of particles is produced. Although the collinear factorization breaks 
down, the perturbation theory does not due to smallness of the coupling at the scale Qs even 
for the IR modes. Failure of the collinear factorization of parton distribution functions has 
been discussed already in the pioneering publications on gluon saturation . It was sug- 
gested that a more general type of factorization, fc-r-factorization, may hold at high energies. 
However, it turned out that although the fc^-factorization is much better approximation of 
the exact formulas than the collinear factorization, it is violated as well in all high energy 
processes save for the single inclusive gluon production in 7* A collisions jl^ . 

As the consequence of breakdown of factorization of the fragmentation functions, 
hadronization pattern changes with parton density in both target and projectile and thus, 
exhibits complicated energy and atomic number dependence. In the present paper we dis- 
cuss one of the possible hadronization channels which explicitly breaks down the collinear 
factorization. From empirical point of view, failure of the collinear factorization is evident 
from observations of strong energy dependence of baryon to meson ratios in pp and dAu 
collisions jl6]. 

Another process in which the collinear factorization breaks down is the relativistic heavy- 



In fact, factorization theorems break down at even higher scale, leading to the so-called extended geometric 
scaling, see Ref. 0, 0, [l3|. 
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ion collisions in which the produced system of color charges probably evolves through the 
dense and hot stage (Quark Gluon Plasma). However, it must be heeded that the collinear 
factorization is violated already at the early stages of decay of the classical fields, preceding 
the formation of the QGP. It has been suggested that parton recombination may be an 



alternative mechanism of hadronization 17, 



18 



19 



2l|. It yields a surprisingly good 



description of experimental data on elliptic flow and might also be relevant for interpretation 



of baryon to meson ratios reported in 



151]. It is based on a simple idea that small-fcr niesons 



are formed by coalescence of two constituent quarks with transverse momenta /ct/2 at the 
same rapidity, while baryons are formed by coalescence of three constituent quarks with 
momenta at the same rapidity. 

Motivated by the phenomenological success of the recombination approach we set to 
investigate the process of gluon recombination in the high gluon density regime. In this 
paper we consider pion production and argue that it is dominated by the recombination of 
two classical fields. Our approach naturally incorporates momentum conservation as well as 
the recombination geometry in the coordinate space (recombining partons must be in the 
same elementary volume of phase-space). In heavy-ion collisions it yields a "cold" nuclear 
matter effect on particle hadronization and may be essential for analysis of the residual 
"hot" nuclear matter effect. 

The recombination process which we discuss in the present paper consists of two stages: 
(i) production of a pair of gluons in J^^ = 0'*'^ and color singlet state; (ii) recombination of 
this pair into a pair of vr-mesons using the anomaly matching mechanism 22|]. The relevant 
diagrams are depicted in Fig. [TJ Shown in Fig. [TJ^a) is emission of two gluons by valence 
quarks belonging to different nucleons in the incoming nucleus, their successive merging and 
production of a gluon pair in the singlet J^*" = O"*"^ state. 

An alternative way to produce a pair of gluons in a singlet J^'" = 0~^~^ state is shown in 
Fig. [l](b). Unlike the diagram Fig. [it^a) it represents a first step in quantum evolution which 
proceeds by emission of a soft gluon. This type of evolution gives rise to the Kharzeev- 



Levin soft Pomeron 



fsee also 



23 



24|). However, this diagram is parametrically small as 



compared to Fig. [I](a). Indeed, the diagram of Fig. [T](a) is of the order of ctg Qt qI ^ ^, where 
Qt and Qp are the parton densities in the target and projectile respectively; Qt ~ A^^^, where 
A is the atomic mass of the target nucleus (if the projectile is a nucleus of atomic mass B, 
then Qp ~ B^^^). Here we used the fact that in a quasi-classical (McLerran-Venugopalan) 
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(c) 



FIG. 1: Examples of diagrams contributing to the production of a pair of gluo ns in J^^ = 0"*"^ 



and color singlet state: (a) quasi-classical case, (b) first radiative correction 



22 



23|, 



24( 1 ■ (c) single 



gluon production followed by the conventional fragmentation [25|]. Horizontal solid lines are the 
valence quarks belonging to a different nucleons. Vertical dashed line describes an instantaneous 
interaction which can happen at different light cone times (for simplicity we show here only one 
possible interaction; see Fig. [2] for other possibilities). 

approximation a1 Qt ^ a1 Qp 1. On the other hand, the diagram Fig. [T](b) is parametrically 
of the order of Qt Qp ^ as- 

At low invariant masses, description of pair production in the singlet J"^*" = 0+"'" channel 
in terms of color degrees of freedom becomes inadequate. Spectral density of the corre- 
sponding correlator (see Sec. IIIip is saturated by colorless excitations the most prominent 
of which are pions. Unlike gluons which contribute to the spectral density at the order 
Og, see ( l36l) . pions contribute at the order see( H2i) . Therefore, diagram Fig. [T](a) is of 

ly l arg er than 
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28| 



the order of a'^ Qt ~ 1/ a1 at low invariant masses. This is parametrica! 
the hadron production via the coUinear fragmentation of a single gluon Q, 
shown in Fig. [I](c). Indeed, the corresponding single gluon production diagram is of the 
order of Qt Qp ~ I /as- On the other hand, diagram Fig. [l](a) represents a higher twist 
effect as compared to the single gluon production and thus has an additional suppression 
factor Ql/k"^ at high transverse momenta (here Qs is associated with the projectile). These 
two hadronization processes become of the same order at transverse momenta of the order 
of k"^ ~ Ql/as- We expect that at lower transverse momenta the recombination mechanism 
discussed in this paper gives the main contribution to the particle hadronization at high en- 
ergies. At RHIC energies, the corresponding kinematic region is about A;^ < 3 GeV for light 
hadrons in Au-Au collisions at midrapidity. It is significantly wider at forward rapidities 
and at higher energies. 
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The paper is organized as follows. In Sec. Ullwe calculate the double inclusive cross section 
for glu on p air production in the singlet J^^ = 0^^ channel in the framework of the dipole 
model 33|] . Since the relevant degrees of freedom at low invariant masses are pions which 
have large inelastic cross section on a nucleon, we are going to neglect the diagram Fig. [2p 
in which the produced pair interacts with the target. Indeed, the survival probability of 
a pion in a heavy nucleus is exponentially suppressed as compared to the one of a color 
dipole which has much smaller characteristic size ~ ^/Qs- This approximation amounts to 
assumption that the intermediate gluons in Fig. [2] are almost on-mass-shell, see (|3]). The 
resulting "wave function" is given by ( 12^ . We then calculate the forward amplitude of each 
color dipole in a nucleus, ( JT9l) and that of the entire projectile ( I2T1) . In doing that we neglect 
correlation between partons belonging to different nucleons since the corresponding dipoles 
have sizes about 1.3 fm (typical nucleon separation). The double inclusive cross section of 
gluon pair production is the convolution of the projectile "wave function" and the forward 
scattering amplitude in the coordinate space and is given by ( !26l) . 

In Sec. mil following in steps of 221], we describe the formalism of anomaly matching 
and calculate the double inclusive cross section for vr-meson production, ( H3|) . Our approach 
is based on the observation that the scale anomaly is closely related to the finite density 
of vacuum fluctuations of quantum gluon fields. These fluctuations are characterized by a 



semi-hard scale Mq ~ 2 - 2.5 GeV |29|, 



30| . The presence of this scale makes reasonable the 



perturbative expansion and allows for the calculation of non-perturbative contributions to 
the spectral density due to the scale anomaly. Our method is similar to the QCD sum rules 
approach. The later is based on the operator product expansion which holds only in the 
presence of a hard scale. 

In Sec. |IV]the quantum evolution effects are taken into account. This makes it possible to 
address energy, rapidity and atomic mass dependence of the pion production cross section. 
The final result is given in Eq. ( !54l) . 

The fragmentation process suggested in this paper has a number of phenomenological 
consequences. Firstly, the energy dependence of the cross section is steeper than in the 
conventional fragmentation mechanism since it requires exchange of an additional Pomeron 
between the projectile and the produced pair. Secondly, as Fig. |4] implies, at higher energies 
heavier hadrons can be produced along with pions. Therefore this hadronization mechanism 
yields a non-trivial energy dependence for different particle species and can be used to 
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FIG. 2: Contributions to production of a gluon pair in J^^ = O"*""*" and color singlet state in a 
quasi-classical approximation. Only diagrams in which the gluon at X2 is emitted from a valence 
quark before the gluon at xi are shown. 

analyze the energy dependence of particle ratios. This issue is addressed in a great detail 
in the forthcoming publication [31]. Other two-gluon channels besides the J"^*" = 0+^ may 
also have interesting phenomenological applications and can be considered along the same 
lines. 



II. TWO-GLUON PRODUCTION IN A QUASI-CLASSICAL APPROXIMATION 



In the light-cone perturbation theory 32j there are eight diagrams contributing to the 
"wave function" of a gluon pair. In Fig. [2]four of them are shown; the other four diagrams can 
be obtained by switching the order of gluon emission from valence quarks. Note that in the 
eikonal approximation, the life-time of a parton fluctuation in the fast nucleus wave function 
t{ ^ fc+/k^ is much larger than the typical time of interaction with the target U ~ Ra, 
where kj^ is the large light-cone momentum of the parton, k its transverse momentum and 
Ra is the nuclear radius. Therefore, in this approximation we can regard the interaction as 
instantaneous. This observation constitutes the basis of the dipole model jssl. It implies 
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that an inclusive production cross section is a convolution of a projectile "wave function" 
with the dipole scattering amplitude in the transverse coordinate space. 

Diagrams in Fig. [2] differ one from another by the structure of their energy denominators 
which we would like to consider now more closely. Each of the eight diagrams in Fig. [2] 
contains three energy denominators. In the diagrams A-C the two leftmost energy denom- 
inators corresponding to emission of gluons at Xi and X2 are different for each diagram 
whereas the rightmost one, corresponding to the 2— >-2 gluon scattering process is the same 
in all cases. The product of the first two energy denominators is given by 

A +A ■ — ^ + — 1 - fia) 

' ki_ + k2- ki_ ki^ + A;2_ ki_ /c2- ' 

B ■ ^ 1 _ 1 1 ^^^^ 

k2- k2- + p'L - v'- k2- k2- - h- - h- ' 

■ k^ ki_+p'L-p'_ ^ k^ A;i_ - /i_ - k- ' 

1 1 _ 1 1 

21- p'i _ p>_ k2- + p'L -v'-~~ h- + k- k2- - h- - h- ' ' 

while the last energy denominator is 

1 1 



(If) 



ki_ + k2- + p'L - pL ki_ + fca- - h~ - h- 

where we used the overall energy conservation condition p"_ + /i_ + k^ = p'_. The product 

of the energy denominators in the diagram D is 

/ 1 1 \ 1 1 1 1 

^21 + ^12 : \ = • (Ig) 

V^i- ^2-7 A;i_ + k2- h- + k- ki^k2^ /i_ + 12- ^ ^' 

Simple calculation shows that sum of all energy denominators from all eight diagrams van- 
ishes. We conclude that in the absence of interactions 

(^^. + + + ^A.) = (2) 

« J=l,2 

as expected. 

As explained in the Introduction, the correct description of the particle pair produced at 
low invariant masses in the scalar J^^ = 0~^~^ state is furnished using the hadronic degrees 
of freedom. By virtue of the color transparency, pion-nucleon inelastic cross section is much 
larger than the dipole-nucleon one. Hence, we are tempted to neglect the contribution 
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of diagram D to inclusive cross section. This imposes a certain constraint on the energy 
denominators A-C since we must ensure that the condition ([2]) holds, for otherwise the 
resulting cross section would not be gauge invariant. In other words, we require that 

E (^A,+^B,+^cJ ^0. (3) 

This approximate equation must hold with the same accuracy as the assumption that pions 
are completely absorbed by the nucleus. It can be guaranteed if the intermediate gluons are 
almost on mass-shell, i. e. ki- + ~ h- + h-- Then, Eqs. ((H) become 

A +A ■ — ^ + — ^ - — — (4a) 

' k2- ki^ + k2- /ci- ki^ + k2- /ci- k2- ' 

„ „ 1 1 1 1 11 

L-21 + C-12 : — ■ — J} + — ■ — Y, ~ - — . (,4cj 

— pL /C2- + P- — P- P- — P~ ki^ -\- p_ — p_ /C2_ 

Obviously, sum of all energy denominators in (jl]) conforms to condition ([3]). 

Since the structure of the energy denominators in cases A-C is the same, we need to 
calculate the "wave function" in only one case, say, A. Introducing momentum q such that 
(see Fig. [2]) 

ki = ^ih + l2)-q, k2 = ^{h + l2) + q, (5) 



we have 



A1A2 



u{pi) ^A,. u{pi-{h + l2)/2 + q) 1 
ViPi - (^1 + ^2)/2 + q)+ ki- 



c^2* u{p2-{h + l2)/2-q) 1 



V(P2 - ill + /2)/2 - q)+ h- 

y \ Ai A2pM!^po- / A3* A4*\ r 

^ + k2- - h- - I2- - 20 ^ " -'"^ ^ " 

((/i + l2)/2 - g)+ ((/i + l2)/2 + g)+ IGvrS ' ^ > 

where V^^^ is the four-gluon vertex: 

(7) 

and (. . .) g means projection onto the singlet state. After projection of the final gluons onto 
the color singlet state, the color factor becomes the same for all six terms contributing to 
the four-gluon vertex. 
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In the eikonal approximation, we can simplify the q ^ qg vertex as follows 33 1 



— 7 ■ ^ — . = - — ~ , loj 

Vp^ ViPi - ^i)+ '^i 

where the bold typeface distinguishes the transverse component of the corresponding four- 
vector. It is convenient to introduce the following notations: (i) fractions z and z' of the 
light-cone momenta of each of the four gluons are given by 

h+ = ^^"^ ^ -q+ = z ih+ + I2+) , k2+ = ^ ^^"^ +q+ = {1- z) + I2+) . (9a) 

/i+ = z' + k+) , I2+ = (1 - z') + I2+) ; (9b) 
(ii) the total and the relative transverse momenta k, k and q are defined as 

k = li + l2 = ki + k2, = z'h-H- z')\i, q = ;zk2- (l-2)ki. (10) 

Note that the invariant mass of the produced pair is 

= (/i + hf = + h-l2+ - 2 li • I2 = \ . (11) 

Z [I z ) 

With the help of these equations, after some simple algebra, we can derive 

1 1 



+ k2- - k- -k^-iO I _ JL _ '2 _ ,n 

zk+ ^ {l-z)k+ z'k+ {l-z')k+ ^ 

k+zil- z] z' il - z') 



(12) 



q2 z'il - z') - K,^z{l - z)-iO 

(We added expression — (ki + k2)^/A;+ + (li + l2)^/A;_|_ = to the denominator of the first 
line) . 

At high energy, the scattering matrix is diagonal with respect to the color dipoles. There- 
fore, in order to take the interaction of the incoming parton system with the nucleus into 
account we need to transform the wave function \I'(/i, I2) into the coordinate representation. 
This is accomplished as follows: 

*te.«J-/gi/(g.^-"-^'--^^*a..w- (13) 

Here ^1 and ^2 are the coordinates of the final gluons with momenta Ii and I2, respectively. 
It is convenient to introduce the transverse coordinates of the intermediate gluons xi and 
X2. To this end, we change the set of integration variables in (fT3|l and from {Ii, I2, q} to 
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a new one {k, q} using ([5]) and (fTOj) . The Jacobian of this transformation is unity. The 
phase factor in f|T3l) becomes 

-^h■^l-^h■^2 = -^^^■{^2-^l)-^ (ki + k^) ■ (/^ + (1 - 6) • (14) 

We can identify xi = X2 = z' + (1 — z') $,2 = x. We observe that the coordinates of 
both intermediate gluons are equal. This is the resuh of the integration over the internal 
momentum q. Indeed, had we Fourier transformed the amplitude with respect to li, I2 and 
q, the subsequent integration over q would have given the delta function (5(x2 — xi). 

Summation over the gluon polarizations can be carried out using the rule {|e^e^*) = 6ij/2. 
Recalling that in the light-cone gauge e ■ 1] = = we derive using ([7]) 

^11 ^ abed ^fj /5 1*^1 ^ J 1^2 ■ e )dcd 

Ai,A2 

= -g'Yl ■ ^^') + ■ e"^*) (k2 ■ e^^*) 5x,x, 

Ai,A2 

= - / (ki ■ k2) irj''^ + r'^f^'] 5,, Sx,x, = -g' (ki ■ k2) 2 AT, Sat <5a3A, . (15) 
We can now write the "wave function" as 

^-(^,x) ^ -.^T.T.2iV.....3A. I ^2 I |^e-~/ §,l\z 

(^k-q)-((l-z)k + q) z'{l-z') 

(zk-q)2((l-^)k + q)2 q2z'(l-z') -k2^(1 -;z) -iO ' ^ ^ 

where ^ = - 

The cross section for the double-inclusive gluon production is given by 

(Ij0^gg X 1 _ _ ^ ^ f tO 1^ f -tO f -lO f -lO f dZ 



d%d%dyd% 8Ar2 

A3A4 

x^^3A4(^^ x)*^«^^*(r7, y) 2 H(x, y) e*ii-(''-y+(i-^') (^-'')) e^^^'^^-y-"' , (17) 

where the coordinates y and 77 in the complex-conjugate amplitude correspond to the coor- 
dinates X and $, in the amplitude; S(x, y) is the rescattering factor, b is impact parameter 
of the projectile with respect to the target and the color traces are taken over each nu- 
cleon (hence the subscripts 1,2). Introducing color gluon notations in the amplitude as 
shown in Fig. [2] and denoting by letters with bars the corresponding gluon colors in the 
complex-conjugated amplitude, we find the overall color factor as 

— — Tr(7~' T ) Tr(X' (^f^^^ jbde _j_ jade J'^'^^^ ^jace jbde _j_ jade jbce^ ^cd 

1 ^aa ^bb ^ab ^ j^^^ab^ ^ab ^ j^^^ab^ ^ 3 CpN^ . (18) 



4iV2 
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We can see that all four non-vanishing gluon indices permutations in the four-gluon vertex 
give the same color factor. 

Next, we turn to calculation of the forward scattering amplitudes for each time and gluon 
ordering. Let us for a moment assume that the intermediate gluons have an arbitrary coor- 
dinates xi and X2 in the amplitude and yi and y2 in the complex conjugate one. Calculation 




AA 







AB 



AC 




FIG. 3: All possible time sequences for interaction of a projectile with the target. The diagrams BA, 
CA and CB can be obtained by complex-conjugation of the diagrams AB, AC and BC respectively. 
Interaction of the final gluon pair is neglected as explained in the text. 



of the scattering amplitudes is similar to the case of a single gluon production [25|. We have 





= e 


i(xi-yi)2Q2 ^_l(x2-y2)2Q2 ^ 




'^BB 


= e~ 


\{^i-yi?Ql ^ g-i(x2-y2)2(; 


n + e-i(x?+yi)Q? + g- 




= 1, 






Sab 


= e~ 


ix2Q2 ^_l(x2_y2)2+Q2 ^ 


ix2Q2 ^_l(x,_y,)2+Q2 ^ 


Sac 


= e~ 


i„2r)2 _i„2r,2 






= e~ 


l„2r)2 lY2n2 
_j_ g 4-'^2^s 





(19a) 
, (19b) 
(19c) 
(19d) 
(19e) 
(19f) 

Taking into account Eqs. Pa|) - fj4c|) we can sum up all the scattering amplitudes to obtain 

^AA + ^BB + ^CC — ^AB — ^BA + ^AC + ^CA — ^BC + ^CB = 
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In the limit of recombining gluons Xi = X2 = x and yi = y2 = y we derive 

S(x,y) = (l + e-3(--y)'«' -e-^^'«' -e-^y'«')' . (21) 

To write the final expression for the double-inclusive gluon production, we introduce an 
auxiliary function F{^,x) such that 

^A3A4(^^ x) = - Ta n 2N, 6ab Sx,x, F{^, x) . (22) 

With this notation we obtain 

"^""^^ - N^Cf [ (fi [ (fx [ dPr, [ fy ^ 



d^hd%dyd% 27r2 " y V J J J z'{l-z') 
xF{^, x) F*{r], y) 2 H(x, y) e^ii-(^-y+(i-^') («-^)) g^i^-Cx-y-^' (4-r,)) _ ^23) 

Eq. (l23l) is a general result which we derived in the eikonal approximation and assuming 
that the recombining gluons are almost on-mass-shell. It can be significantly simplified if we 
note that the recombining gluons must be close in rapidity which implies that the light-cone 
momentum fractions carried by gluons are typically equal z,z' ~ 1/2. Assuming that this 
configuration is dominant, we derive in Appendix A the following expression for the "wave 
function" 

^^3A4(^^ ^) ^ 2iV, 5ab 5x,x, ^ In M - — j (24) 
and the cross-section 

dagg 2ag NcCp f d^x f d^y 



dHidH2dyd% vr^ (I2 

X [1 + iy\h - 111 K^{-iy\\2 - U|)]e*('i+i^)-(--y) 2 s(x, y) . (25) 



Using (12511 we can extract the inclusive cross section for production of a pair of gluons 
with invariant mass M. For this objective, it is convenient to consider the inclusive cross 
section in terms of momenta k = I2 -|- li and /t = (I2 — li)/2. Then integrating over the 
directions of k and recalling that by (ITT]) = M'^/A yields 

dagg _ N^Cp f d^x f d^y ^ik.(x-y) 



dM^d'^kdyd% tt^ J x^ J y^ 

X [l-ixMKi{ixM)][l + iyMKi{-iyM)]2E{x,y). (26) 
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Integrating^ over k we get 

= — " ; ' / —\l-txMKi(txM)\ 8 l-e"4^^. . 27) 

The integrand of (!27|) can be expressed in terms of Bessel functions exploiting the relation 

|l - ixMKi{ixM)f = (l + I MxYi{Mx)y + y (Afx)^ /^^(Ma;) . (28) 

In the limit of large invariant masses M ^ Qs the dominant contribution to the integral 
in fl271) comes from dipoles of size 1/M <^ s <^ 1/Qs- Expanding both the McDonald 
function and the dipole scattering amplitude yields 



M:^Qs. (29) 



dM^ d% SvrS 

In the opposite limit of small invariant masses the largest (logarithmic) contribution stems 
from sizes 1 /Q^ <^ x ^ 1 /M in which case we estimate 



da„„ GAa^N^Cp 1 



772' M<^Qs. (30) 



dM^ d% 7r6 M2 

This behavior has important phenomenological consequences which will be elucidated in the 
upcoming publication |31 |. 

Finally, the total inclusive cross section is determined from (123!) by first integrating over 
n which yields the delta function 5(^ — 77). Then using fICll) and flA8P we derive 



d'^b Svr^ Jq X 3tc 

Numerical calculation of the ratio of (|27|1 to (|3T1) is exhibited in Fig. HI 



III. ANOMALY MATCHING 



Now as we derived the cross section for the double-inclusive gluon production in the scalar 
color singlet channel J^*" = O"*""*", we can use the anomaly matching procedure to derive the 
double inclusive pion production at low invariant masses. This approach has been discussed 
in details in |22j]. Here we give a brief review. 



^ Convergence of the integral in (pS)) is discussed in Appendix C. 
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FIG. 4: Ratio R = -^^^j/^^ (l^) function of invariant mass M. Values of parameters 

for solid line: Qs = 1 GeV, dashed line: Qs = 2 GeV. In both cases e/x = 1 GeV. 

In the product of the "wave function" and its complex conjugate, there appears a loop 
formed by the produced gluons. The contribution of this loop is proportional to the corre- 
lator'^ 

(0|T{^;:(x)^,^(0)}|0), (32) 
where 6'^'^ is the energy-momentum tensor. In the chiral limit its trace acquires a finite value 



^' = -^2^'"''^ (33) 



due to the scale anomaly of QCD. The correlator (15^ can written in the spectral represen- 
tation 

m') = ^ J d'xe^'- (0|T{^^;^(x)C(0)}|0) = j da' , (34) 

where the spectral density 

Peik') = Yl I '^^{2nf5^'\pn-k)\{nr^m\ (35) 

k'^ = is the invariant mass of the produced system. In the lowest order of perturbation 
theory the spectral density is given by 

Pr(M^)^(|l)^H^MV (36) 



^ All coordinate and momenta notations in this section are independent from the notations in other sections 
unless otherwise specified. 
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Although the scale invariance of QCD (in the chiral limit) is broken down by quantum 
fluctuations, there remains a residual symmetry which manifests itself in an inflnite tower 



of equations, known as the low energy theorems 29|], relating various Green's functions 
involving operator 0'^{x). The flrst in this tower of equations relates the Green's function of 
the flrst and second order as follows 

n(0) = -4(0|^|:|0) = -16evac. (37) 

Spectral density (!35l) represents sum over all physical states, perturbative (high M^) and 
non-perturbative (low M^). We can see using (!36|) in (IMI) that the perturbative contribution 
to n(0) is divergent. Therefore, it must be subtracted to satisfy the theorem fl37|) : 



[pf^^(a2)-pr(^')] = -16 6vac. (38) 



(t2 



Thus, the vacuum expectation value of 6'|^ measures the energy density of non-perturbative 
fluctuations of vacuum. 

If the invariant mass of the produced gluons is small, we can no longer describe the 
produced particles in terms of the color degrees of freedom. Rather it is appropriate to 
express O'^ directly in terms of hadrons. This can be done using the effective chiral Lagrangian 



40| 



£ = ^ Tr d,U d^U^ + ^ml fl Tr(t/ + f/t) , (39) 

where U = e^*'^/-^'", vr = rr'^t"' and t"" are the SU{2) generators. The trace of the energy- 
momentum tensor for this Lagrangian is 

01^ = -d^TT" 9^7r° + 2 m^TT'^TT'' + . . . . (40) 

In the chiral limit we have 

{n^n-\9^,\0) = M\ (41) 
leading to the following non-perturbative contribution to the spectral density 

PViM') = ^,M\ (42) 

where is assumed to be less than a certain cutoff Mq at which the perturbation theory 
becomes applicable. This cutoff is related to the vacuum energy density Cvac and is estimated 



to be rather large Mq ~ 2 — 2.5 GeV 



29 
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The main idea of anomaly matching is that while at large M ^ the physical spectral density 
coincides with the perturbative formula fl36|) . at low it is specified by f l42|) . According 
to (l38l) gluons do not contribute to the spectral density at low M^. The two formulas 
must coincide at the scale Mq. Therefore, in order to calculate production of two pions at 
M < Mo we need to replace the perturbative contribution to the spectral density (l36l) — 
calculated (indirectly) in the previous section — by the two-pion contribution (H2|) . This 
substitution is equivalent to calculating the diagram in Fig. [51 Therefore, the cross section 





1 












K 



FIG. 5: Production of a pion pair in J^^ = 0^"*" channel with invariant mass M < Mq. 
for the double-inclusive pion production becomes (using notation of the previous section) 

dM'^d?kdy(Ph ~ 7r5 62 J J ^ 

X [l-ixMKi{ixM)][l + iyMKi{-iyM)]2E{x,y), (43) 

where we employed (!36l) and (H2|) in (126|) . It is remarkable that this cross section is of the 
order a^, two remaining factors of Og are inherent to the gluon distribution functions of each 
projectile nucleon. 



IV. INCLUDING QUANTUM EVOLUTION 

R 

Here we are going to include the small-x nonlinear quantum evolution of 51 into the 

n 

cross section from Eq. fl26l) . Since the evolution equations in [5| are written for the forward 
amplitude of a quark dipole on a nucleus, we have to first generalize Eq. (12^ to the case 
of gluon pair production in two dipoles-nucleus scattering. Of course, such a model of 
a nucleon is a rough approximation. However, in the gluon saturation regime details of 
nucleon structure play a little role. Therefore, our results below may still serve as a good 
approximation to a more accurate treatment of the nucleon jl^. The generalization of (126]) 
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to dipole-nucleus scattering is easily done by including emissions of the s-channel gluon 
in Fig. [2] by the quark and anti-quark in the incoming dipoles. Denote the transverse 
coordinates of the quark and anti-quark in the incoming dipoles by "zq and "zi, where the 
superscript a = 1,2 labels different nucleons. Then, instead of ( l26i) we have 

_^^3^_n 2 ^ _ 8aJ iV^ f<fxf A .k.(x-y) V V 2^fx v"z- «z ) 



1 - i|x - "zi|Miri(i|x - "zi|M) 



where "zqi = "zq — ''zi and 



l + i|y-%-|Miri(-i|y-%-|M) 



-i(x--z,)2Q2 



(44) 



(45) 



The inclusion of quantum corrections in the leading logarithmic approximation (resum- 
ming powers of UsV) in the large- Ai'c limit is done along the lines of [l^ using Mueller's 



dipole model formalism 



33|. Since we assume that the produced gluons are at the same 



rapidity, the prescription for inclusion of quantum evolution is identical to the single gluon 
production case. We first define the quantity ni("zo, "^zi; "wq, "^wi; F — y), which has the 
meaning of the number of dipoles with transverse coordinates "^wo, "^wi at rapidity y gener- 
ated by the evolution from the original dipole ""Tiq^Tii having rapidity Y . It obeys the dipole 



equivalent of the BFKL evolution equation 



33 



38| 



gni(°Zo,"zi;"wo,"wi;|/) 
dy 



with the initial condition 



2 7r2 



d''Z2 



-01 



^20 ^21 



niC'zo,%;''wo,"wi;y) 



+ni(%, "zi; ''wo, "wi; y) - ni("zo, "zi; '^Wq, "wi; y) 



(46) 



ni(%, '^wo, "^wi; y = 0) = 5("zo - ''wq) 6{^z, - ^w,] 



(47) 



If the target nucleus has rapidity 0, the incoming dipole has rapidity Y, and the produced 
gluons have rapidity y, the inclusion of small-x evolution in the rapidity interval Y — y is 
accomplished by replacing the cross section from fl44p by [l^ 



d a, 



gg 



dW^ d^kdyd% 



'^Zoi,^Zoi) 



f H d'^wo rf'^wi ni(%, '^zi; ^Wo, ''w,; Y - y) 



d a, 



gg 



dW^ d^kdyd% 



;^woi,Voi). (4J 
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Eq. (HHj) neglects correlations between the original dipoles as explained in Sec. HI The 



evolution in each of the original dipoles is linear as was originally shown in [ij] : the Pomeron 
splittings cancel in the rapidity interval between y and Y in compliance with the AGK cutting 
rules [39 1 . 

Inclusion of evolution in the interval between and y is accomplished by replacing the 
Mueller-Glauber rescattering exponents according to the following rule 14 1 



-i(xo-Xi)2Q2 



l-iV(xo,xi,F) 



(49) 



where A^(xo, xi, Y) is the forward amphtude for a quark dipole xq, xi scattering on a target 
with rapidity interval Y between the dipole and the target. It obeys the following evolution 
equation 



9A^(xo,xi,y) asNc f ,2 



dY 



2 7r2 



d X2 



^01 



-^21 



[iV(xo, X2, Y) + iV(x2, xi, Y) - iV(xo, xi, Y) 



-7V(xo,x2,r)iv(x2,xi,r)] 



with the initial condition 



iV(xo,xi,y = 0) = l-e-^^^o-''^)' 



(50) 



(51) 



Performing the substitution from fH9|) in psl) yields 



S(x, y; "z,, ^z; Y) = (iV(x, Y) + N{y, Y) - iVC^z,, '^z,, Y) - iV(x, y, Y)f . (52) 



With the definition of Eqs. (1521) we write the following answer for the double inclusive 
gluon production cross section in the scalar J^*" = O"*"^ channel including small-x evolution 
effects 

2 



da. 



99 



( zoi, ZOly 



8aJ N,Cf 
~ M4 



7r 



j n d-^wo rf^wi mC^Zo, '^zi; "^wo, "^wi; F - y) 



dM2 d'^kdyd'^h 

^ fd^ f ^^.Hx-y) ^ 2H(x/w.;y/w,;y) 

J J y' 

X [l-i|x-''wi|Miri(i|x-"wi|M)] [l+i|y-''wj|MJs:i(-i|y-"wj|M)] . (53) 
Similarly, the double inclusive pion production cross section is given using ( l43l) by 



a=l 

1 



dM2 d^kdyd^h 



;^zoi,^zoi) 



96 a? 1 



552 



/" n ^"wo t/"wi mC^Zo, "zi; '^wo, '^wi; F - y) 

^ a=l 
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X [l-z|x-"wi|Miri(i|x-"wi|M)] [l + i|y-"wj-|Miri(-i|y-"wj|M)] . (54) 



This is the central result of our paper. We are going to use it for the phenomenological 



analysis of the RHIC data in the forthcoming publication 



3l|. 
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APPENDIX A: DERIVATION OF THE INCLUSIVE CROSS SECTION IN THE 

z,z' = 1/2 APPROXIMATION 

Using expressions ( |T6l) . (I22i) the function F(^,x) becomes 

d'k f <fK _,,.^_,k.x f dh , (|k-q)-(lk + q) 



' J (27r)2 y (27r)2 J IQti^ ( i k - q)2 (i k + q)2 q2 - ,^2 _ " ^ ' 

Integration over k can be done using the following formula (see Appendix B for the deriva- 
tion) 

/■ (ik-q) ■ (ik+q) ..^ Iq-x , , , 

/ (0 ^2 )\ = — ^ sin 2q-x . A2 

J (27r)2 (ik - q)2 (ik + q)2 vrq^x^ 

The subsequent integration over At is a Fourier transformation of a two-dimensional Feynman 

propagator and can be expressed in terms of the Hankel function 

"'"e- ^, \_ .„ -^gftoO. (A3) 



(27r)2 q2 _ ^2 _ ^0 4 

where q= |q|, ^ = |^|. With the aid of ( 1A2I) and (1A3I) we obtain 



• poo 

F(|,x) = / dqH^^\qO M2qx) . (A4) 

This integral can be taken by analytically continuing the Hankel function to imaginary values 
of q with the help of Eq. 9.6.4 of Ref. ^ 



^o(^r) = -yi^f (r), (A5) 
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and using Eq. 6.576.3 of Ref. [35|. The result is 

1 



To calculate the cross section (1231) we use the following integral 



(A6) 



d2^e-^5«(i2-ii)in(i_ 



4x" 



27r 



d^^Jo 7,ah-h\] In 



16 vr 



— \l - ix\h - 111 Ki{ix\h - li 

I2 - li)'' 



(A7) 



Inserting flA6l) and flA7p into fl23|) and using K*{iz) = Ki{—iz) we derive fl25l) . 
Let us also note for future reference the following integral 



In 1 



2x 



27l^X^ 



+ 2'k'^x^ - 



+ / <f«ln ( 1 



(A8) 



APPENDIX B: DERIVATION OF (fA2]) 



Consider an auxiliary function G'(x, q) defined as following 



G(x,q) 



One can readily verify that 



d^k (|k-q)-(|k + q) 
(27r)2(lk-q)2(ik + q)2 



-i k x 



ik-x 



[ik-q)2(ik + q)2 



(Bl) 



-U^x + ^q -V.-^q G(x,q). (B2) 



Thus, the problem is reduced to evaluation of (^(x, q). Using the Feynman's trick we write 



G(x,q) 



da 



d'^k 
(2^ 



jkx 



Hik-q)2 + (l-a)(lk + q)2}2 



Introducing a new vector k' = ^k + (1 — 2a) q we come by 



G(x,q) = 4^'rfae'(^-2")*'i"y" 



d^ 



-2jk'-x 



1 



X 



^^g2(l-2a)^q.x_ 



{k'^ + 4a(l-«)q2}2 
Kl{Axq^^a{l - a)) . 



(B3) 



Now, the integral (IBip is dominated by two IR logarithmic singularities at k = ±2q 
which tie in with emission of soft gluons ki = and k2 = by the projectile. In (1B3p 



21 



they correspond to the values a = 0, 1 of the Feynman parameter a. Keeping only the 
logarithmically divergent terms we derive 



(B4) 



T^q Jo 

Changing the integration variable a = jS"^ yields 

G'(x,q) = — cos(2q-x) lim / d(3 (3^ Ki{Axq(3) = — cos(2q-x) ln(l/gx) . (B5) 

Finally, using flB5l) in flB2p we arrive at flA2p . 

An important remark is in order here. Eq. (1A2I) is not valid in two limits: (i) x — > 0, q 
fixed and (ii) q ^ 0, x fixed, though it holds in the limit qx —>■ 0. The reason is that the 
integral on the left-hand-side of (lA2p is symmetric with respect to transformation x ^ q. 
This can be verified using twice the formula 

2^5 (™) 

to transform the integrand into the coordinate space. This symmetry property is violated 
in the above mentioned limiting cases (i) and (ii). The values of the integral in these cases 
are 



k 

k2 



and [251 



d^k 



~tkx 



87rln(l/a;/x) 



(ik-q).(|k + q) 
(ik-q)2ak + q)2 



]7T ln(g//i) , 



where is an IR cutoff. 



APPENDIX C: CONVERGENCE OF THE INTEGRAL IN (l26l) 

Consider the following integral 
I{M,k) = -^ J^J ^e^''-(^-y)[l-2xMiri(^xM)] [l + iyMK,{-iyM)]E{^,y 

The dipole scattering amplitude S(x, y) < 1 by unitarity. Therefore, 

(fx ^ r°° dx 

——- e'^-"" \1 - ixM KMxM)] = / — Mkx) \l - ixM KAixM)] 
(27r) x^ ^ Jq X -' 



/(M, k) < 
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where x = |x| and k = |k| as usual. Since (see (1A7P ) 

y ^jQ{kx) [1 - ixMKiiixM)] = In ( — - ij 



we derive 



/(M, k) < 



In 



M2 

"12' 



- 1 



(CI) 



(C2) 



This expression is finite apart from the logarithmic divergence at k = M corresponding to 
li ■ I2 = 0. Integration over M in (126|) is thus bounded by (see (lASP ) 



dM 
M3 



I{M,k) < 



dM 
M3 



In^ 



M2 
1^ 



+ 



In' 



M2 

A;2 



1 



57r2 
6fc2 



(C3) 
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